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Abstract 

We present two new results on relativistic hydrodynamics with anomalies and 
' external electromagnetic fields, "Chiral MagnetoHydroDynamics" (CMHD). First, 

. we study CMHD in four dimensions at second order in the derivative expansion 

assuming the conformal/Weyl invariance. We classify all possible independent con- 
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formal second order viscous corrections to the energy-momentum tensor and to the 
U{1) current in the presence of external electric and/or magnetic fields, and identify 
eighteen terms that originate from the triangle anomaly. We then propose and mo- 
tivate the following guiding principle to constrain the CMHD: the anomaly-induced 
^ ■ terms that are even under the time reversal invariance should not contribute to the 

^ , local entropy production rate. This allows us to fix thirteen out of the eighteen 



transport coefficients that enter the second order formulation of CMHD. We also 
relate one of our second order transport coefficients to the chiral shear waves. Our 
second subject is hydrodynamics with (A^ + l)-gon anomaly in an arbitrary 2N 
dimensions. The effects from the (A^ + l)-gon anomaly appear in hydrodynamics 
at {N — l)'th order in the derivative expansion, and we identify precisely such 
corrections to the U{1) current. The time reversal constraint is powerful enough to 
allow us to find the analytic expressions for all transport coefficients. We confirm 
the validity of our results (and of the proposed guiding principle) by an explicit 
fiuid/gravity computation within the AdS/CFT correspondence. 
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1 Introduction and Summary 



Quantum anomalies are among the most beautiful, subtle, and important phenomena in 
quantum field theory. Recently it has become clear that anomalies play a very important 
role also in the macroscopic dynamics of relativistic fiuids. Much of this progress is 
motivated by the possibility to observe the anomalous "chiral magnetic" currents in non- 
Abelian quark-gluon plasma produced at Relativistic Heavy Ion Collider and the Large 
Hadron Collider. However, anomalous currents in hydrodynamics can be expected to 
have many other applications: for example, the time-reversal- invariant, non-dissipative 
currents in strongly correlated systems that we will discuss in this paper present clear 
interest for quantum computing. 

The Chiral Magnetic Effect (CME) is the anomaly-induced phenomenon of electric 
charge separation along the axis of the applied magnetic field in the presence of fiuctuating 
topological charge [H 121 |3l HJ [5] . The CME in QCD coupled to electromagnetism assumes 
a chirality asymmetry between left- and right-handed quarks, parameterized by a chiral 
chemical potential fiA- Such an asymmetry can arise if there is an asymmetry between the 
topology-changing transitions early in a heavy ion collision. Closely related phenomena 
have been discussed in the physics of neutrinos ^ , primordial electroweak plasma [7] and 
quantum wires |8j. 

While the original derivation used the weak coupling methods (see also Refs.[9l 
for recent discussions on this), the origin of the effect is essentially topological and so 
the CME is not renormalized even at strong coupling, as was shown by the holographic 
methods [IH [121 [131 |ll|15l HSl [IT] . The evidence for the CME has been found in lattice 
QCD coupled to electromagnetism, both within the quenched approximation [HI [191 120] 
and with light domain wall fermions [21j. Unlike the baryon chemical potential, the chiral 
chemical potential fi^ does not present a "sign problem" which opens a possibility for 
lattice computations at finite [4] . A direct lattice study of the chiral magnetic current 
as a function of fiA was performed very recently [22j; it confirms the expected dependence 
of CME on the chiral chemical potential and the magnetic field. 

Recently, STAR [231 [21] and PHENIX [23 [26] Collaborations at Relativistic Heavy 
Ion Collider at BNL reported experimental observation of charge asymmetry fiuctuations. 
While the interpretation of the observed effect is still under intense discussion, the fluc- 
tuations in charge asymmetry have been predicted [1] to occur in heavy ion collisions due 
to the CME. Very recently, STAR reported the expected [3l [1] disappearance of the 
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effect at low collision energies where the energy density of created matter is smaller and 
likely below the critical one needed for the restoration of chiral symmetry. The ALICE 
Collaboration at the Large Hadron Collider at CERN has just reported [28] the observa- 
tion of charge asymmetry fluctuation signaling the persistence of the effect at very high 
energy densities. Additional future tests of CME include the positive correlation between 
the electric and baryon charge asymmetries [29] . See also Ref . [30] . 

Since in the strong coupling regime the plasma represents a fluid (for reviews, see [3T1 
[32]). it is of great interest to study the effects of anomalies in relativistic hydrodynamics. 
A purely hydro dynamical derivation of the anomaly effects at first order in derivative 
expansion was given by Son and Surowka [33|, motivated by earlier results in AdS/CFT 
correspondence [3ll [351 El] which found, among others, chiral vortical effect. It has been 
generalized to anomalous superfiuids [371 ES [39] and non-abelian symmetry [lOj [H] . It 
has been found that the CME current persists in hydrodynamics [12] and is transferred 
by the sound-like gapless excitation - "the chiral magnetic wave" [43l [44], see also [45] 
for an earlier study of collective excitations in anomalous hydrodynamics. Another line 
of development has been along the viewpoint of the effective field theory [16] . 

The idea of Son and Surowka [33] was to consider the local entropy production rate 
d^s^ and to impose on it the positivity constraint following from the second law of ther- 
modynamics. The contributions from the anomaly to the entropy production were shown 
to be locally unbounded in either sign so that unless their coefficients identically van- 
ished, they could potentially violate the second law of thermodynamics. These arguments 
lead to a set of algebro-differential equations for the transport coefficients related to the 
anomaly; in many cases they can be solved. 

The present work continues the investigation of anomalies in relativistic hydrodynam- 
ics. We also consider the entropy production as an important constraint. However, we 
propose a different guiding principle: instead of requiring the positivity of the total en- 
tropy production rate, we argue that the anomaly-induced viscous corrections should 
not contribute to the entropy production at all, due to the time-reversal invariance of 
the anomalous transport coefficients. The time-reversal T invariance provides a unique 
criterion that can be used to establish the nature of currents. For example, the "usual" 
electric conductivity a is T-odd, as can be easily inferred from the Ohm's law J' = aE^: 
the electric field is T-even, whereas the electric current J* is T-odd. On the other hand, 
the (anomalous) quantum Hall conductance is a T-even quantity, as it is associated with 
a T-odd magnetic field. The physical meaning of T invariance of transport coefficients is 
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quite simple: T-odd conductivities describe dissipative currents, whereas T-even conduc- 
tivities describe non- dissipative currents. The anomaly-induced currents are protected by 
topology and are thus of non-disipative nature; as such, they do not contribute to the 
entropy production. We will discuss this in more detail in section 12. 4( we will verify that 
our guiding principle leads to non-trivial relations among the transport coefficients that 
are obeyed by the available results of the explicit holographic computations in section [231 

Let us add that enforcing a positivity constraint on the entropy production is not easy 
when one considers a second (or higher) order in the derivative expansion in hydrody- 
namics. Indeed, one first has to consider only a subspace of possible configurations with 
vanishing previous order contributions to meaningfully discuss constraints at the second 
(or higher) order. This typically results in very few useful constraints on the second or 
higher order transport coefficients. However, the T-invariance guiding principle that we 
propose here provides a stronger constraint on the anomaly-induced terms, and in some 
cases allows to evaluate them. Note that the first order hydrodynamics is known to have 
problems with causality and to be numerically unstable, so from a practical point of view 
the second order formulation of relativistic hydrodynamics with anomalies and external 
electric/magnetic fields is highly desirable; it is the topic of the present paper. 

The paper is organized as follows. In the first part (section [2]), we consider the 
anomaly-induced viscous corrections at second order in the derivative expansion in four 
space-time dimensions. As the number of possible independent terms increases drastically 
at the second and higher orders, we assume the underlying conformal symmetry to con- 
strain the problem. We first classify all possible second order viscous corrections to the 
energy-momentum tensor and the U{1) current, including also the non-anomalous terms, 
in the presence of external electric/magnetic fields. Let us mention that Refs. [Ml 135] pre- 
viously classified possible second order conformal viscous terms without including external 
electric/magnetic fields. Then, by considering discrete symmetries of charge conjugation 
and parity, we select eighteen terms with transport coefficients that are necessarily linear 
in the anomaly coefficient. These transport coefficients are time-reversal T even, and we 
demand that they do not contribute to the entropy production. This enables us to fix 
thirteen out of the eighteen anomaly-induced transport coefficients. 

Our results for the second order transport coefficients related to the triangle anomaly 
are new and can be expected to be universal. Four transport coefficients of interest for us 
were previously computed in the holographic AdS/CFT approach using the fiuid/gravity 
correspondence [3ll [35] , which makes it possible to check and confirm some of our results. 
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This test is quite non-trivial, and we consider it as an important evidence for the vahdity 
of our proposed guiding principle and of the relations that we derive from it. However 
we also identify many other second order transport coefficients coming from the anomaly 
that have not been computed in the AdS/CFT setup or in any other framework, and a 
more thorough fluid/gravity computations for them are certainly very desirable. 

We also explain the phenomenon of the chiral shear wave, that is, a helicity-dependent 
shear mode dispersion relation due to triangle anomaly, in terms of one of the second order 
anomalous transport coefficients in the energy-momentum tensor. Chiral shear wave 
was first observed in Refs. [171 SHI SS] via linearized hydrodynamic analysis in AdS / CFT 
correspondence. However its relation to viscous transport coefficients was unclear and we 
close this remaining gap in the present paper. 

In the second part of this work (section [3]), we consider the anomaly effects in hy- 
drodynamics in a higher 2N dimensional spacetime, where the theory has an underlying 
(A^ + l)-gon anomaly. By using the charge conjugation and parity symmetries, we show 
that the effects of the anomaly first appear at (A — l)'th order in the derivative expansion 
in the f/(l)/entropy currents, and we identify precisely A such terms. Although these 
terms are of a very high order in derivatives, the time reversal invariance still dictates 
that they should not contribute to the entropy production, and this principle provides us 
a sufficient set of constraints to determine these terms completely and derive for them 
analytic expressions. We then confirm our results via the fiuid/gravity correspondence in 
an AdS/CFT setup, corroborating our guiding principle. 

2 Second order relativistic conformal hydrodynamics 
with triangle anomaly 

2.1 A primer on the confer mal/Weyl covariant formalism 

In this subsection we review the formalism of conformal hydrodynamics basing mainly on 
Refs. [501 EI] ; this will also allow us to introduce the notation. Conformal hydrodynamics 
by definition is covariant under Weyl transformations. 




(2.1) 
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where is an arbitrary scalar function on the spacetime. In our convention, a Weyl 
covariant tensor with Weyl weight w transforms as 

, (2.2) 

so the metric tensor has w = —2. Note that the upper indexed inverse metric g'^'^ has 
w = +2 instead. In short, the constraints from conformal symmetry on hydrodynamics are 
simple: i) the energy- momentum tensor and the local symmetry currents that constitute 
the basic elements of hydrodynamics should be Weyl covariant; ii) the energy- momentum 
tensor should be traceless up to a local Weyl anomaly (which in most cases is a higher order 
effect in the derivative expansion scheme). One can easily derive Ward-type identities for 
the assumed Weyl invariance, for example as in Ref.[50], and obtain the transformation 
properties of the energy-momentum tensor and the currents. The energy-momentum 
tensor T^'^ which is traceless = has a Weyl weight w = +6 and a current has 
w = +4. In hydrodynamics, one is dealing with locally varying thermodynamic variables 
such as temperature T, pressure p, chemical potentials /x, local velocity u^^, etc, and writes 
the energy-momentum tensor and symmetry currents in terms of these variables - the 
resulting expressions are the so-called constitutive relations. Imposing Weyl covariance 
on these expressions is a powerful constraint that reduces much arbitrariness one might 
have without the conformal symmetry. In writing constitutive relations, one typically 
invokes the derivative expansion scheme; the zeroth order expressions in the Landau 
frame are 

T^"" = {e + p)uf'u'' + pg^"" = Apu^'u'' + pg"" + ■ ■ ■ , (2.3) 
f = nu^' + --- , (2.4) 

where we have used the tracelessness condition e = 3p in the first line; n is the charge 
density. This gives us the Weyl weights of p {w = +4), n {w = -|-3), and {w = +1). 
The temperature T naturally has w = +1 and the Gibbs-Duhem relation e + p = Ts + fin 
implies that the entropy density s has w = +3 and the chemical potential fi has w = +1. 
The zeroth order expression for the entropy current = su^ + ■ ■ ■ implies that s'^ has 

W = +4:. 

A useful mnemonic used in Ref . [50] is that the Weyl weight of a Weyl covariant tensor 
is given by 

w = [mass dimension] -|- of upper indices] — of lower indices] . (2.5) 
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For example, an external gauge potential that couples to has to = 1 — 1 = 0, which 
is also confirmed by considering the Weyl invariance of the action 

Sr^ I d^x^A^f . (2.6) 

Let us point out yet another useful nmemonic derived from holography. A tensor in 
the 4 dimensional conformal field theory appears in the holographic 5 dimensions as a 
component of a 5 dimensional bulk field. The radial wavefunction of that particular 
component near the AdS boundary has a form r""" when the AdS metric is written as 
^+r'^{dx^y. One can identify w in the exponent as the Weyl weight of the corresponding 
tensor. For example, a bulk gauge field contains a current and an external gauge 
potential as 

4^^^' = ^^^ + ^. ; (2-7) 

this expression conforms to the fact that has w = +2 (recall that = g^uj") and A^ 
has w = 0. The 4 dimensional metric gf^,y is present in the bulk metric as 

9lT = r'g,u , (2.8) 

which is also consistent with w = —2. This mnemonic is natural because the holographic 
versions of Weyl transformations are in fact the coordinate reparameterizations r — )■ e^'^r. 

For higher derivative viscous terms, the construction of Weyl covariant expressions out 
of derivatives becomes non-trivial, simply because the "ordinary" covariant derivatives 
are not Weyl covariant. The Weyl covariant formalism from Ref. |5T] that we will 
use is based on the idea of introducing Weyl covariant derivatives which contain 
suitable Weyl connections to make them covariant under the Weyl transformations. This 
is in close analogy with the usual electromagnetic/gravitational covariant derivatives. In 
the case of electromagnetism, the connection introduced is a new dynamical degree of 
freedom one has to add to the theory, but in our case of hydrodynamics we do not want 
to introduce new degrees of freedom. Therefore, one needs to construct a Weyl connection 
that transforms properly out of the already existing hydrodynamic variables; in our case it 
is the velocity field m^. The situation is quite similar to the metric Christoffel connection 
that is constructed out of the metric itself. We refer the reader for details to Ref. [51], and 
will outline only the main features of this approach. 

The Weyl covariant derivatives share almost all common properties with the co- 
variant derivatives, including the chain rule. The V^T^Z has the same Weyl weight w 
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of the original T^.. '. so that "D^ has w = 0. Note that has w = +2. The important 
properties of for the purpose of hydrodynamics are 

u^VX = , ^^ = , V^g,^ = , V^e'"'^^ = , (2.9) 

where e""'^'^ is the totally anti-symmetric covariant tensor which is easily found to have 
w = +4. To illustrate how acts, acting on a scalar / of Weyl weight w is 

V,f = V,f + wW, , (2.10) 

where W/^ is an analog of electromagnetic connection, but constructed out of as 

= n^V.n^ - ^^n^ , (2.11) 

and w is acting as a charge in electromagnetic analogy. Note that is first order in 
the derivative. For tensors with indices, Vfj_ involves not only wWfj,, but in addition more 
connections acting on tensor indices like metric connections. We will rarely need this 
detail, but an important fact is that these connections are all linear in so that 
increases the number of derivatives by precisely one. As expected, commutators of I^^ 
bring us several kinds of Weyl covariant curvature tensors. The simplest one is 

[V^,V,]f = wW^J , (2.12) 

where W^^jy = d^yVy — dyW^ is a new Weyl covariant tensor with weight tli = 0. Another 
example is 

[V^,V,]V^ = wW^,V^-'R^JVp , (2.13) 

which includes a Weyl covariant cousin of the Riemann tensor; however its symmetry 
properties are slightly different [5T] . 

The main point is that the basic hydrodynamic equations are in fact Weyl covariant 
- that is, one can show that 

v^T^"" = v;r^^' , v^f = v^f , (2.i4) 

precisely when T^'^ {j^) has w = +6 {w = +4) and is traceless = 0. Therefore the 
Weyl covariant formalism can present a useful framework for the studies of conformal 
hydrodynamics. 
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2.2 First order conformal hydrodynamics with anomaly 

We will begin by applying the Weyl formalism to the first order hydrodynamics with 
triangle anomaly. The strategy is closely parallel to the one presented in Ref.[53]: one 
imposes the positivity condition on the local entropy production, which turns out powerful 
enough to allow the determination of the transport coefficients. Compared to the original 
study in Ref. |33] of the general non-conformal case, we will see that the derivation in the 
conformal case has some subtle differences, although the results will be identical. This 
first order case is a warm-up exercise prior to a more elaborate study of the second order 
hydrodynamics in the next subsection. For simplicity we consider a single U{1) current. 

The basic Weyl covariant equations of hydrodynamics are 

= ^e^''''^F,,F^p = KE''B, , (2.15) 

o 

where F^y = d^A^ — d^A^ = V^Ay — V^A^ is the field strength tensor of external gauge 
potential A^, which has a Weyl weight w = and we define 

E^' = F^''uy , = ]^e^^'^^UyF^p , (2.16) 

both of which have w = +3. They are also transverse, E'^Ufj, = B^u^ = 0. One can check 
that (12.151) including the anomaly is consistent with Weyl weights. We then write the 
constitutive relations of T'^" and j'^ in terms of the local thermodynamic variables, 

Tf"" = Apu'^u" + pg>"' + T^,] + T^^^ + ■ ■ ■ = Apu^u" + pg>"' + Tf"" , 
f = nu'' + u^^^ + u^^) + ■ ■ ■ = nu^' + ly^ , (2.17) 

where the subscripts in the viscous terms denote the number of derivatives each term 
contains. In the Landau frame defined as 

T''''uy = 3pu'' , jf'u^ = -n , (2.18) 

the viscous terms must be strictly transverse order by order, 

^rn> = ^(V/. = . (2.19) 

Following Ref. [323 consider F^^ as of first order in derivative, which is a weak field 
limit. 
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Among the thermodynamic variables {p^T, ^i^u^) with u^u^ = —1 or their combina- 
tions, there are five independent ones that can be chosen to describe a local property of 
the plasma; for convenience, we choose two scalars (T, /i = ^) and u'^. The number of 
equations in fl2.15p is also five and the time evolution of the system is well-defined. In 
writing down possible viscous terms containing derivatives, one can use the basic equa- 
tions of motion in f l2.15p to replace certain first order derivative terms with higher order 
derivatives, effectively removing them at a given fixed order in derivatives. Explicitly, the 
first equation in fl2.15p with fl2.17p gives 

4 {V^p) uV + Vp + P^r^'^ = nE'' + . (2.20) 

Contracting with li^, one obtains 

u^V^p = ^ {E^v^ + u.V.T^'') , (2.21) 

and inserting this into the original equation, one gets 

V^p = nEf" + F^"z/« - ^u" {E^u^ + u^V^r'"') - V^t^" . (2.22) 



Similarly, from the second equation in ( I2.15p . one obtains 



u^'V^n = -V^pf" + K Ef'B^ . (2.23) 

Now, any scalar thermodynamic variable should be a function of two variables we choose, 
i.e. (T, fi); note that T has the Weyl weight w = +1 and fi has w = 0. This means that 
a scalar / with a weight w can always be written as 

f = T-f = T^m ; (2.24) 

we will use this quite often. Writing p = T^p and n = T^n, one sees that eqns. (I2.22p . 
f l2.23p can be used to remove P^T and u^V^Ji in favor of E^ and Af^^V^fi up to higher 
derivatives, where A^,^ = {g^y + u^Uy) is the projection operator to the space transverse 
to u^. Therefore, one can always ignore them in writing down possible viscous terms at 
a given order in the derivative expansion. This means that the only derivative term of a 
scalar quantity one needs to consider in constructing the viscous terms is A^j^yV^Ji] indeed, 
any scalar is a function of {T,p) only, and the only other possible derivatives are V^Uy. 
This is a big simplification. It will also appear useful that i?^ can be replaced by 

Efj^ = —V^p + higher derivatives . (2.25) 
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Since the number of equations in (12.151) is five and we have used them to remove five com- 
ponents V^T and u^V^li, no further reduction is possible from the equations of motion. 

Following Ref. [33], one considers the local entropy production rate V^s^ = V^s'^ 
where the equality is valid since has w = +4, and imposes the positivity condition on 
all possible configurations, which will lead to a few differential equations for the transport 
coefficients. It appears possible to solve them using conformal symmetry. As usual, one 
starts from 

= (4 {V^p) u^u" + + V^T^ + ^1 {{V^n) + V^p^) 

= -u^V^{3p) + fiu^V^n + u.V^T^" + ^V^y^ 

= -u''V^e + ^u^'V^n-(V^u,)r^' + fiV^iy^ , (2.26) 

using (12. 9p . e = 3p and u^rf^^ = 0. From the thermodynamic relations de = Tds + /udn, 
e + p = ip = Ts + fin, and the Weyl weights of each quantity, one can check that 

V^e = TV^s + fiV^n ; (2.27) 

and using this relation in the above formula after some algebra gives 

TV^ {su'' - fMu'') = - (V^u,) t''" - {TV^Ji - E;) v'^-C fiE'^B^ . (2.28) 

One proceeds by writing the entropy current s'^ in the derivative expansion as 

s^' = su" - Jiu" + sfi) + sf2) + ■ ■ • , (2.29) 

so that the total entropy production rate we want to keep positive definite is 

ra^s^ = -a^.r'^'^ - (TP^A - E^) u'' - k fiE'^B^ + TV^ (sf^) + sf^) + ■■■)> , (2.30) 

where we define 

o-^,u = ^{T^f^Uu + T^uUfj.) , ^^fiu = '^(Pf.Uu -T^uUfj.) , ^^^lUy = + u^^ . (2.31) 

Note that a^i, and uj^j_i, have w = —1 and are already tansverse and traceless due to (12.91) . 
The constraint (12.301) is the main starting point in considering the entropy production; it 
can be rewritten in a simple form within the conformal formalism. The task is to classify 
^ j^fi^ order by order in derivatives and to find constraints on the transport 

coefficients associated to them by imposing positivity on the local entropy current (I2.30p . 
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At first order in derivatives, independent available tensors are 

A^^'V.fi , V^u, = a^, + u^, , ; B'' , (2.32) 
it is also useful to consider a transverse pseudo-vector of w = +2, 

_ l^/..a/3^^^^^ ^ le^-/3^^v„w;3 . (2.33) 

For the first order energy-momentum tensor r^'^^ which must be transverse and traceless, 
there is only one possibility as can be deduced from f l2.32p : cr^'^. Taking into account the 
positivity of entropy production f l2.30p it must assume the classic form 

rir^ = -2r^a'^'' (r/ > 0) , (2.34) 

where the shear viscosity r] has a weight w = +3. 

For z/^^^ and s^^^, one needs to construct the transverse vectors and there are four 
of them: A^'^V^Ji, E^', B'', and u^". The positivity of fl2:30|l allows to fix easily the 
dependence on the first two: 

u^,) = -a{TA^''V,Ji-E^) + iuj^ + iBB'' , (2.35) 
sfi) = Du^' + DbB'' , (2.36) 

with a positive conductivity a of weight w = +1\ the remaining transport coefficients ^ 
{w = +2), C,B iyj = +1), D {w = +2), and Db {w = +1) are to be determined. We 
will rigorously show in the next subsection using discrete symmetries that they should be 
proportional to the anomaly coefficient k, as they originate from the anomaly, but for now 
we will simply let them be general possible terms as in Ref . [33j ; we then insert the above 
into fl^TSUj) . 

To proceed, we need the following two identities to be valid at all orders: 

V^u^ = , V^B^^ = -2E^io^ . (2.37) 
Proof : From (12.91) and fl2.3ip one has 

= e>''"'^u^,u^fs + ^e>''"'^u,[V^,V^]ufs . (2.38) 
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Because u^i, is transverse, that is, in the local rest frame of = the only non-vanishing 
components are ojij = 1,2,3), the first term clearly vanishes. For the second piece, 
use the relation fl2.13p 

[D^,V^]up = -W^^Uf,-n^Jus , (2.39) 
and a symmetry property j51] 

T^,J + T^^M=0 , (2.40) 

so that 

2V,u^ = -le>^--Pu,n^Jus = -\e^'^^^u, {u^J + TZ^^^J) us = . (2.41) 
For the second identity, one has 

= e'^^-^c^^.F^^ + e^'^^^M, (V^F„^) . (2.42) 

The second piece vanishes via Bianchi identity, while for the first piece it is most convenient 
to go to the local rest frame where 

21)^5^ = 2e^^'=Vii^fco • (2.43) 
Noting that in this frame we have 

E, = F,ou'' , 0;^= = ^e^o^^iioo;,, = -^wV^'^V, , (2.44) 
one finally concludes that 

2P^S^ = -AEku'' = -AE^oo,, . (QED) (2.45) 



Inserting (^^3^, flOGl) to flOOj) and using (^M), one arrives at 

TP^s^ = 27]a^,a^'' + a{TA>^''V,fi-E^){TA^^V"fi-E^) 

- {TV^fi - E>^) i^u^ + CbB^) - K ^lE^B^ 
+ T{{V^D)io^-2DBE^uj^ + {V^DB)B^) 

- a,,T^^.-{TV,fi~E^)u^,.+TV^s^,. + --- , (2.46) 
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where the first hne is manifestly positive definite and the last line is of higher order. In 
considering arbitrary configurations, it is important to remember that is not indepen- 
dent, but is equal to ^l^fiP up to higher derivative corrections. Writing p = T^p{fi), this 
implies is given in terms of T>^T and T>^fi. However, some other terms in the above, 
V^D and V^Db, are also expressed in terms of P^T and T>nfi as we write D = T^D and 
Db = TDb, and therefore, E^ is not completely independent of V^D and V^Db as one 
considers arbitrary configurations at the first order. The easiest systematic way to deal 
with this subtle difference from the non-conformal case is to simply replace i?^ with ^l^^p 
by (12221): 

E^' = -V^'p - -F^"z/„ + — (E'^Ua + u^V^r'"') + -I)„r"^ , (2.47) 
n n on n 

at a given derivative order in considering arbitrary configurations. One then finds after 
some algebra 

TV^s^' = 2r]a^,a^'' + a{TA^^''V,fi-E^){TA^^V"fi-E^) 
+ (^-nV^f, + TV,D +{^- V^^ uj^ 

+ l-nBV^fL + TV^Db +(—- /€- ) V^p] B'' 



n n 



— OuvT, 



{TV,U-E^)v^^.+TV^s^^^.+ 



+ -{-F>^"u^ + V^T"^^){{^-2TDB)uj^ + i^B-^if^)B,) ; (2.48) 
n 

the last line is a remnant from replacing i?^ by ^l^^iP using f l2.22p . which is in fact relevant 
when we discuss the second order viscous terms in the next subsection. 

Because u'^ and B^^ can take arbitrary values and also can change their signs, the sec- 
ond and third lines can easily overcome the first unless the coefficients vanish identically. 
This leads to two differential equations 

Ct 2TD \ 
n ^)^^^ " ° ' ^^-^^^ 

-T^bV^U + TV^Db +( — - K^] V^p = . (2.50) 

\n n J 

Let us now substitute 

p = T'p{fi) , {D,0 = T\D,0 , iDB,CB) = nDB,^B) , n = T'n , (2.51) 
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upon which the above becomes 

T'(-^+D' + ^[^-2DB)p'^V^fi + T'(2D + ^[^-2DB)'ip^V,T = 0, 

where prime denotes As T and ft are independent, the coefficients in front of V^T 
and T>^fi must vanish separately, so that one in fact gets four equations to solve. 

We now need to prove the following relation: 

p' = n . (2.52) 

Proof : Let us start from the basic thermodynamic relations e + p = 4p = Ts + fin and 
de = 3dp = Tds + /idn, which give 

Ap = s + pn , (2.53) 
3 (AT^pdT + T^p'dp) = T {3T^sdT + T^'dp) + Tp (3r%dT + T^n'dp) . 

As T and p are independent variables, the coefficients of dT and dp on both sides in 
the second relation should agree separately. The equation from dT is precisely the ffist 
relation, while from dp one obtains 

3p' = s' + pn' = Ap' - (pn)' + pn' = Ap - n , (2.54) 

where we use the first relation in the second equality. One gets p' = n from the final 
expression. (QED) 

Using fl2.52p . the four equations simplify as 

D' + 2Db = , D'b-kP = , (2.55) 

— — fl — — fl — 

i = 2DB--D , ^j, = ^j,- Db , (2.56) 
2p Ap 

from which it is easy to integrate D and as 

D = \np'' , Db = \kP'' , (2.57) 
and the second line finally gives 

One can check that our results agree precisely with the general non-conformal results in 
Ref . [33] upon using the conformal relation e + p = Ap. 



*We are neglecting possible integration constants ^AT , which might be related to gravitational anoma- 
lies [10]. 
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2.3 Second order conformal hydrodynamics with anomaly 



In this subsection, we address our main objective of studying the second order viscous 
corrections to the energy- momentum tensor and the current, with particular attention to 
the anomaly-induced effects. Our starting point is Eq.f l2.48p from the previous subsection: 



+ (-T^V^fi + TV^D + - V^p^ 
- - (TV^Ji - E^) 1/(^2) + r^^A.sfa) + " " " 

+ l(_FMV,(^)+I),r('^5')((e-2rDBK + (eB-«:/x)5,,) . (2.59) 

The first three lines are already taken care of in the previous subsection; for example 
the second and the third lines simply vanish when we use our solutions for the first order 
transport coefficients ^, ^s, D, Db- The last two lines are what is important in this subsec- 
tion. We are interested in possible second order viscous corrections r^'^'^, z/|^j, s^2)i ^^^1 by 
considering the entropy production we would like to constrain the transport coefficients 
associated with them as much as possible. Although we will classify all possible inde- 
pendent second order viscous corrections, we only focus on the transport coefficients that 
necessarily arise from the triangle anomaly in our consideration of entropy production. 
In other words, we are going to specify to a subclass of possible terms with transport co- 
efficients that are linear in the anomaly constant k. Due to a selection rule from discrete 
symmetries of charge conjugation C and parity P that we will discuss shortly, there is 
no interference in f l2.59p between the anomaly-induced terms that we focus on and other 
"non-anomalous" terms, so that one can meaningfully separate them in the consideration 
of entropy production. We leave for the future the task of fully exploring the constraints 
from entropy production on all possible "non- anomalous" transport coefficients that we 
list. 

Let us begin with z^^^-j and s^^) ^^'^ write down all possible independent Weyl covari- 
ant transverse vectors that include two derivatives. As discussed before, hydrodynamic 
equations of motion can be used to remove covariant derivatives of scalars such as P^T 
in favor of and A^'^Vj^p,. Other possible tensors of use are T>f^Ui, and B^. It is quite a 
tedious job to list all possible combinations one can construct and to identify independent 
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components, so we simply present the resulting fifteen terms: 

af'^E^ , a^^B^ , u^'^E^ , uj^"" ^u^V^Ef" , (2.60) 
e^-'^^'^u^E^VpfL , e^'^'^u^B^V^fi , e^'^'^u^E^Bp , e^'^^^u^V^E/s , e'^^'^^ u,V ^B p . 

The five terms in the first line were found before in Refs. |34[ [35] , and the rest are new. 
There are a few details involved in showing that these are indeed all possible independent 
terms. For example, u^V^B^ is not included in the above because it is related to the term 
e^'^^^uJDaEij by Bianchi identity. Another possibility is to use the second order Weyl 
curvature tensors to construct terms such as W^'^Uy and A^^7?.j,q,u", but one can show 
that 

V'V^u, = (a^, + uj^,) = -W^^u" - n^.u" , (2.61) 
and moreover, using the relation 

one can easily check that W^^Wy and A^^TZuaU" are already included in the above list. 
We leave other checks to readers. 

As we are interested in terms that are necessarily linear in the anomaly coefficient k, we 
need a systematic method to identify such terms. We will use two discrete symmetries; 
charge conjugation C and space parity in the local rest frame P. Consider the basic 
hydrodynamic equation 

V^f = t,E^B, . (2.63) 

Under {C,P), the spatial component of is (—1,-1). Similarly, E'^ has (—1,-1) and 
B^^ has (— 1,+1) 0. The covariant derivative has (+1,-1). This tells us that n has 
(C.P) = (—1,-1). From the constitutive relation = nu^ + ■ ■ ■, the charge density 
n and the chemical potential Ji have {C,P) = (— 1,+1). The combination (k/I) then 
has {C,P) = (+1,-1). Naturally, the entropy current has {C,P) = (+1,-1). The 
usefulness of these discrete charges is exposed by the following observation: when we write 
linear combination of the terms in the above list, the transport coefficient 
in front of each term should have a well-defined (C, P) that is easily derived by comparing 

'''We define our P such that the gauge potential transforms as — ^ and -> —A^. Alter- 
natively, one can introduce an additional overall negative sign as in the case of axial symmetry in QCD. 
It is a matter of interchanging P o CP, and as long as one adheres to a chosen definition, it would not 
make a difference for our purpose. 
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the (C, P) of iy^2) ^{2) "with that of each term in the hst. These transport coefficients are 
conformal scalars of some weight w, so that they can be generically written as T'^f{p,, k). 
Since they have definite {C,P) and /t, /t/i) all have different (C, P)'s, there are only 



four possibilities of the form of /(//, k); 






(CP) = (+1,+1) : 


■ f{l^,i^) 




(CP) = (-1,+1) : 


■ f{l^,i^) 




(CP) = (-1,-1) : 


■ f{l^,i^) 


= Kgip,"^, K^) 


{C,P) = (+1,-1) 


■ f{l^,i^) 





Therefore one can systematically select the terms that necessarily come from the anomaly 
by choosing only the terms whose transport coefficients have (C, P) = (±1,-1). In 
retrospect, the first order transport coefficients $,,^b, D, Db have precisely such {C,P) 
charges. 

With the help of this criterion, we find that five terms in the above list can be identified 
as originating from the anomaly: 

a'^'^u, , a^'^B, , u^^B, , e^^^'^u^E^Vf.ft , e^'^'^u^V^E^ , (2.64) 

and we introduce ten transport coefficients associated with them as 

<2)anomaiy = D^a^^uJ, + D^a'^" B, + D^u^" B, + D.e^'"'^ u^E^V f^fi + D^e'^^'^^u^V^Ep . 

(2.65) 

We perform a similar procedure for rj^^) which has a somewhat more complicated 
structure. Refs. [SH |35] listed all possible independent terms without including external 
electric/magnetic fields, and it is not difficult to extend their results including exter- 
nal electromagnetic fields. Defining the projection operator to transverse, traceless and 
symmetric components as in Ref. [35], 

= ^ (a^A^ + AX - ^A-'^A,,) , (2.66) 
the independent possible terms are 

W^^pB'^V^Ji , W^'pE'^E^ , W^pE'^B^ , n^^P^P^ , W^pe^'^'^a^^u.Es , liy^^'^a^^^u.Bs , 
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where Cp^pcxp is the conformal Weyl tensor of the background metric. Using the (C, P) 
symmetries, one can pick up eight terms that should be hnear in the anomaly k, and we 
introduce transport coefficients for them as 

(2.67) 

The eighteen transport coefficients C,i,Di,Xj {i = l,---,5, j = l,---,8) in fl2.65p and 
( I2.67P are the most general second order transport coefficients of a conformal plasma that 
originate from the underlying triangle anomaly. 



2.4 Time reversal invariance, anomaly and entropy production 

Let us now motivate the main guiding principle that we propose and are going to use 
extensively throughout this paper - namely, that the anomaly-induced T-even terms 
should not contribute to the entropy production. To illustrate the significance of discrete 
symmetries, let us consider first the behavior of contributions to the entropy under the 
spatial parity. The anomalous contributions to the entropy production are special in that 
they change sign under spatial parity transformation P. Suppose there were a contribution 
to the entropy production from the anomalous terms we identify; then in the parity-flipped 
mirror world, this contribution would become negative. Thinking of entropy production 
as originating from some dissipative work, this is very unnatural. This consideration 
gives us the ffist hint that the anomalous terms should not contribute to the entropy 
production. 

The vanishing of the entropy production from the anomaly-induced terms has a simple 
physical meaning - the corresponding anomalous currents are non- dissipative. This rather 
unusual property originates from the fact that the anomalous current is associated with 
the zero fermion modes, and the number of these zero modes is related to the topology 
of gauge fields by the Atiyah-Singer index theorem. Since the topology of gauge fields 
is determined at the boundary of the fluid, the processes in the bulk cannot lead to 
the dissipation of anomalous currents. This consideration can be made more formal by 
considering yet another discrete symmetry of the transport coefficients - time reversal T. 
The "usual" electric conductivity a is a T-odd quantity, as can be easily seen from Ohm's 
law J* = aE^: the electric field is T-even, whereas the electric current J* is T-odd. On the 
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other hand, let us consider the quantum Hall effect as an example of anomalous current in 
(2 + 1) dimensions. The quantum Hall conductance is a T-even quantity, as it is associated 
with a T-odd magnetic field. The corresponding Hall current is non-dissipative, and the 
conductance of the integer quantum Hall effect is given by the Chern numbers of vector 
bundles associated with the energy bands of the Hamiltonian operator [52]. In physical 
terms, the dissipative transport coefficients are described in terms of the response of the 
states near the Fermi energy, whereas the non-dissipative ones involve all of the states 
below the Fermi energy. The anomalous chiral magnetic current can be thought of as a 
quantum phenomenon that involves the entire Dirac sea [53] (reflecting Gribov's view of 
"anomalies as a manifestation of the high momentum collective motion in the vacuum" 
[5ll|55]), and it is thus natural to expect that it is of non-dissipative, reversible nature. 
Indeed, the chiral magnetic conductivity cr^ [56] defined by J = a^B is a manifestly T- 
even quantity as it relates magnetic field and electric current both of which are T-odd. 
We note that this feature of anomalous currents makes them potentially important in 
various applications that include quantum computing, see e.g. [57] . 

We thus conjecture that the terms originating from the anomaly, i.e. the terms that 
are linear in k, do not contribute to the net entropy production at all orders. The first 
order result in the previous subsection obeys this principle - the first order contribution 
coming from anomaly vanishes identically. The validity of this claim in the present second 
order will be evidenced shortly by a non-trivial test of our results against the existing 
holographic computations of some of our transport coefficients. 

2.5 Constraints from time reversal invariance 

Let us now impose the constraints of time reversal invariance (no entropy production 
from the anomaly) on the transport coefficients. We are interested in only the last two 
lines in fl2.59p because the previous lines are already taken care of at pevious orders. 
The last line is an important remnant from the first order computation, and we can insert 
u^^^ = —a {T ts.^^DyJi — E^) and r^'^'^ = —2ria^'^ to get the terms that are linear in anomaly 
coefficient. After some manipulations that use the local rest frame expressions 

F'^Bi = , F'^uj^ = e^'"'^u^B^up = u^''B^ , (2.68) 



19 



the last line in fl2.59p can be written as 

- (e - 2TDb) {TA^^V^fi - E^) u^'^B, - -V, {r]an ((^ - 2TDb) uj, + {^b - B,) . 
n n 

(2.69) 

The most complicated part is to compute to gather the independent compo- 

nents. There are a few non-trivial steps that we have to describe: 

(e^-'^^u^E^Vpfi) = e'^-'-^u^.E^VpfL + e^^^^^w, {V,E^) V^fi , (2.70) 

where we used [V^,Vi3]fi = because /2 has w = 0. In the first term, both u^i, and E^ 
are transverse so that 'Di^fi should necessarily be u^V^Ji which can be removed by using 
equations of motion as we pointed out before; therefore, only the last piece survives. 
Another term is 

V^^e^'^'Pu^V^Ep) = e^''''f'u,,Vo,Ep + ^e^''''('u,[V^,V^]Ep 

= e^'^f^u.^V^Ep + ^e^-'^'^u, {W,^Ep - TZ^JE^) , (2.71) 

and using symmetry of TZ^ap + "^/^lAil] ~ 0, the last term is equal to zero. Also, from 
W^Q, = -^[l^fi^l^alp and Vfj^p ^ nE^ up to higher order derivatives, one can easily show 
that the final result is 

(e^"'''^u,V^Ep) = e^^'^^uj^.V^Ep + -^^^""^u, {V^E^) Ep + higher derivatives . 

(2.72) 

Using the above, the second order entropy production coming from the last two lines 
in fl2.59p becomes 

+\^Il^^'pB^VPji + KW^'pE^B^ + X^W^'pUJ^E^^ 
- [TA^jy^Ji - E^) {iia^-uj, + i2<y^''B, + i^u^'^B, + i^t^^'^^u^V^E^ 

+ {V^Di) e^^'-^u,EJ)pli + {V^D^) e^^^f'u.V^E}^ 
^t{d^ {{V^^o^^) + a^^'V^u,) + D2 {{V^a^^") B, + a^^'V^B,) 
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+D5 (e^^'^'^^co^.Vo^E^ + ^e^^'^'^u, {V^E^) Ep) ) 



+ £ - 2TDb) {TA^^V-fi - E^) u^^-B, - ^V^ {r]an ((^ - 2TDb) + (^b - /^Z^) B,) 



Note that the terms with A3, Ay and ^4 disappear identically, so these transport coefficients 
are simply unconstrained by our method. As before, one should replace E^ by jiT>^p in the 
above, and impose the condition that the total coefficient of each independent component 
should vanish. We list each component and the equation imposed by the condition of zero 
entropy production as follows: 

• From the component a^'^Uv, one obtains 



The transport coefficients are all conformal scalars of some weight, and one can 
generally write them as / = T^f{fi). Inserting this form and noting that P^T 
and T>f^fi are independent, the coefficients in front of these should vanish separately, 
so that the above equation in fact provides two equations. Writing A2 = T'^)^2, 
Ag = TXg, ^1 = T^i, Di = TDi, 1] = T^f] according to their conformal weights, and 
using (12.521) p' = n, one gets two equations as 



(2.73) 



A2'C'/x/U 









(2.74) 
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(2.75) 







where prime as usual denotes 



• From {V^a^^) cOu, one gets 




n 



(2.76) 



• From {V^a^'^) B^, one obtains 




(2.77) 



where D2 = D2 due to its zero conformal weight. 
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From a^'^Bu, one has 

n \ n J n 

which leads to two equations upon writing A5 = TAs, Ae = Ae, ^2 = 

-A5- A6 + 5^--(eB-/t/u)?7' = , (2.78) 

-4A6- + 4e2---(ei?-'^p)^ = . (2.79) 

From [Vfj^uj^'^) B„, one simply gets 

TD3 = . (2.80) 

From u^^By, one has 

{tV^Jx - (^-^3 + ^ (e - 2TDb)) + TV^D^ = , (2.81) 
and using the fact that Z^s = from above, one simply gets 

e3 = -(e-2TDB) , (2.82) 
n 



or equivalently 



1^ = -(1-2Db] , (2.83) 



n 

writing 1^3 = 1^3, = Ta according to the conformal weights. 
From a^yV^u^ , one arrives at 



or upon writing Ai = T^A 



- Ai + TDi = , (2.84) 
Ai = Di . (2.85) 



From a^iyl^^B", it leads to 

- A4 + TD2 = , equivalently X4 = D2 , (2.86) 

where A4 = TA4. 

oj^uV^B" gives one the same equation as in fl2.80p . 

TDa = . (2.87) 
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From e^^^^oj^yVaEp, one simply concludes that 

TDs = . (2.88) 



• What remains can be grouped into two components. One is proportional to 

e^^^'^Ui^VapVpfi and the other is proportional to e^'^^'^UyVa {^^Vj^pj . Now, the latter 
can be expanded as 

e'^^^^w.P, i-Vfip) = -^e'^^'^^u^V^nV^p + ^e^-Z^^^y^ ^ (2.39) 

and the second piece is completely independent of e^'^'^^ UyV ap'D pji, so that one can 
treat the above two components e'^'^^^f^ uJD ^pV j^Ji and e'^'^°'^u,jVa (^^Vf^p^ indepen- 
dently; each component should then vanish separately. The first component gives 

-V^D^ = , (2.90) 

n 

and since Di = D^, it leads to 

= . (2.91) 
The coefficient of the second component leads to 

- ^5 {tV^Ji - ^V^p^ + TV^D, + TD^V^fi + ^D.V^p = , (2.92) 

and using (12.881) = 0, this leads to two equations 

^4 = , 6 = , (2.93) 

where ^5 = ^5. This completes the list of all contraints derived from the requirement 
of no entropy production. 

Let us now solve these equations in a more explicit form. Recall that the first order 
coefficients ^, ^b, Db are already given in the previous subsection in (12.571) and (I2.58p . 
Equations (12T6D . (I2T7D . (ESQ]), (ESH]), (12:851) . (IMj) . (EHH]), and (12:931) trivially give 
solutions for Z)i, D2, -D3, ^3, Ai, A4, 1)5, 1)4, and ^5 as 

= Ai = ^ il-2DB) , D2 = X4 = — (Cb -Kfi) , 6 = - (^-2Db) , 

^3 = ^4 = ^5 = ^5 = . (2.94) 
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From fl2.74p . (12.751) and using the solution for Di, one can solve 

A. + 6 ^ (f(«-2Z)«))V(|-|:)(e-2Z)«) . (2.95) 

As -6 = -^(e-2Ds) , (2.96) 

which fixes A2 + ^i and As — ■^i leaving one unknown such as A2 — i^i. Finally, the equations 
fl2T8|) and fl2T9D with the solution for D2 give 

= (f(e.-K,))' , (2.97) 

A6-e2 = ---(^B-^fi) , (2.98) 

leaving A5 — ^2 undertermined. In summary, we have determined thirteen out of eighteen 
anomalous second order transport coefficients, leaving five parameters A3, A2 — ^1, A5 — ^2, 
Ay, and ^4 unfixed by our method. 

Refs. [M| [35] computed four transport coefficients of interest for us, (Ai^2,35^i) in the 
framework of AdS/CFT correspondence, explicitly corresponding to the = 4 super 
Yang-Mills theory with U{1) R-symmetry at strong coupling. One can test some of our 
relations against these computations; specifically, one can test two relations that we have 
derived, 

A, = H5(4_26,) , 

n ^ ' 

X,H-e. ^ (f(f-2B.))V(|-f:)(f-2«.) . (2.99) 



which are quite non-trivial. Let us first summarize the results of Refs. [34 [ [35], especially 
following notations in Ref.[31]. The 5D action is 



Y 1 MNPQR 

{lQnG5)C5D = R + 12 - -FmnF^^ -—j=^^AMFNpFQn , (2.100) 

4 lzv3 \'~9 

where G5 is related to the gauge theory by = y^. The charged black hole solution is 



m Q"^ 



ds'' = -r^f{r)df + 2dtdr + r^ ^ (^dx') , /(r) = 1 - — + 



1=1,2,3 



^4 3^6 ' 



A = -^dt , (2.101) 
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The temperature T and the chemical potential are given in terms of parameters (m, Q) 

by 

An rfj 

where the horizon rn is the largest solution of firn) = 0. These can be solved explicitly 
as 



(2.103) 



in terms of which the relevant quantities are given by 
_ rh _ ffrU, _ f% 



Sv^TrGsm ' ^ IGv^vrGs ' ^ 24v^7rG'5m ' 

It is satisfying to see that the two relations f l2.99p as well as ^' = fl are obeyed by the 
above holographic results; this is a rather non-trivial test of our guiding principle. It 
would be interesting to perform a full-fledged fluid/gravity correspondence computation 
for other transport coefficients we identify and to check other relations too. 



2.6 Chiral shear wave 

In this subsection, we discuss one physics phenomenon that is related to the second order 
viscous corrections from triangle anomaly addressed in the previous subsections: the 
chiral shear wave. The chiral shear wave is a modification of the transverse shear mode 
dispersion relation 

uj^-i^^k'^±iCk^ + ■■■ , (2.105) 

t + p 

where the leading k"^ piece is as usual, and the term is the first effect originated from 
anomaly. If the system is parity invariant without triangle anomaly, the dispersion relation 
should be invariant under parity transformation A; — )■ — fc, so that only even powers of k 
should have appeared in the dispersion relation. Any odd powers of k are effects from 
triangle anomaly, and the k^ term in the above is the first of them. The ± sign in front 
of it depends on the helicity of the shear modes which will become clear shortly. 
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The chiral shear wave was first observed in Refs. [171 B8| B9] via AdS/CFT correspon- 
dence of = 4 super Yang-Mills with U{1) R-symmetry at strong coupling. The 5D 
action and the charged black hole solution which serves as a background are precisely 
same as in the last subsection, f l2.100p and f l2.10ip . The resulting shear mode dispersion 
relation was |17] 

u;^-iiLe±i e + . (2.106) 

We would like to understand the origin of the piece in terms of our second order viscous 
corrections from triangle anomaly to the energy- momentum tensor and the U{1) current. 
One can easily check that the first order corrections do not induce chiral shear wave, and 
it has been expected that it should be related to higher order corrections §. We will see 
that it actually comes from second order corrections we identify. 

To derive dispersion relations of linearized fluctuations, one starts with a static back- 
ground fluid of temperature T and chemical potential fi, and considers small fluctuations 
of hydrodynamic variables; in our case 6T, 6fi, and 6u'^, and keep only terms that are 
linear in their amplitudes when one writes ST^"^ and 6j'^. For our purpose we don't have 
external electric/magnetic fields, so that the hydrodynamic equations of motion 

d^ST^^'' = , d^Sf = , (2.107) 

will give us spacetime propagation of these linearized fluctuations. In the frequency- 
momentum space, we can read off dispersion relations. 

One first picks up a definite frequency-momentum (u,k = kx^), or equivalently every 
fluctuating mode is assumed to have a common phase factor e~*'^*+*^^'^ . Due to a residual 
S0{2) symmetry in the transverse (x^, x^) space, fluctuating modes are classified by their 
helicities under 5*0(2) rotation, and different helicity modes do not mix. We are interested 
in helicity ±1 transverse shear fluctuations, 

Su±i = (6u^ ± iSu^) . (2.108) 



Note that they are the only modes with this helicity because other possible modes such 
as d2^3{ST,6p.) = are simply absent because 82,3 = by our assumption of momentum 
direction. Therefore, one only needs to consider fluctuations 5u±i. Because we know 
that k^ term in the dispersion relation is coming from anomaly and the usual k"^ term 
from the first order correction, we can neglect non-anomalous second order corrections 



^At least to the authors, it was first pomted out by D.T.Son. Ref.[551 also mentioned it recently. 
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to the energy-momentum tensor and the U{1) current for our purpose. Out of the eight 
anomalous second order corrections to the energy-momentum tensor, one easily sees that 
only the term Xili^pV^ uj^ gives non-zero contributions at our linearized level, so that one 
needs to compute 

5T^'^ = ApSuV + Apu^'Su" - 2r/5a'^^ + Ai^ (n^^P^a;^) , (2.109) 

whose non-vanishing components after some algebra are 

(5r°' = Ap5u^ + ^kuSu^ , 
4 

5T°3 = ApSu^-—ku6u^ , 
4 

5T'^ = -i^kSu^ + ^k^Su^ , 

5T^^ = -i^kSu^ - ^k^Su^ . (2.110) 
Then, the energy-momentum conservation gives 

aoT°2 + diT^^ = -iu (^Ap6u^ + ^ku6u^^ + ik (^-ir]k5u^ + ^k^6u^^ = , 

-103 , amis _ • ,//i„r„,3 -^1 ,f„,2\ , / „-„LX„.3 -^1 l2 e .2 



doT''^ + diT'^ = -iu I ApSu'' - -^ku6u^ \ + ik I -ir]kSu^ - -^k^5u' 1 = 



(2.111) 



-iu(ApTi—ku\+ik\-i7]kTi—k'^\\5u±i=Q , (2.112) 



and in terms of 5u±i = [Su^ ± i6u^), they become 

4 



(2.113) 



which is solved up to order k^ as 



V 1.2 _ . Al ,3 



u ^ -i-^k' T i—k-^ + ■ ■ ■ ■ (2.114) 
4p 16p 

We see that the k^ term indeed originates from one of our second order transport coeffi- 
cients Ai. 

Let us confirm this in the AdS/CFT computation: Ai and p as given by the AdS/CFT 
computations are given in the previous subsection, and recalling Ai = T^Ai, p = T'^p, one 
can check that 

_ A = ^1 _ 4/^^ _ ^h/^^ (2 115) 

IQp IQT^p 24y3T2m2 2AV3m^ ' ^ ' ' 

which agrees with the independent result fl2.106p from the linearized analysis in Ref . |47] . 
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3 Hydrodynamics with anomaly in higher dimensions 



3.1 Entropy current method 

In this second part of the paper, we will relax the constraint of conformal symmetry, 
but instead of staying in four spacetime dimensions, we will explore higher even 2N 
dimensional spacetime with (A^ + l)-gon anomaly among multiple U{1) symmetries. We 
will follow the entropy production constraint similar to the one used in Ref. p3] and in 
the first part of this work. The effects from the anomaly first appear at (A^ — l)'th order 
in derivative expansion, and we will identify N number of terms in the U{1) currents that 
originate from anomaly. We will see that requiring positivity of entropy production fixes 
all the transport coefficients corresponding to them. 

For a concise presentation, we will show details of our derivation for the case of a 
single U{1) symmetry; the results for the multiple U{1) case follow from a straightforward 
generalization of it. The starting point is again the basic hydrodynamic equations 

V -i^ — ^ ^tiivi---mvN p ...p _ DM ('^^^(^\ 



where we define 



DM _ J_ /ii^cn/3r--aiv-i^iv-i,, p ^ ...p „ (^^^7^ 



The meaning of the subscript (0, — 1) in B^^ will become clear in a moment. One 
then invokes derivative expansion in writing down constitutive relations for T^^" and j'^ 
in terms of thermodynamic variables of plasma in the Landau frame as 

f = nu^ + u^" , (3.118) 
where the viscous corrections t'^'^ and u'^ are transverse by the definition of Landau frame, 

Uf.T^"' = u^v^ = . (3.119) 
We assume that the basic thermodynamic relations hold: 

e + p = Ts + , de = Tds + jjidn , dp = sdT + nd^ , (3.120) 
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where s is the entropy density and fi is the chemical potential. As is by now a standard 
procedure, one considers 



= -u^V^e + fiuf'V^n - {e + p - fxn) V/,m^ + u^V^" + fiV^w 
= -TV^ {su^ - fiun - (V^ii.) T^" - T (V^/i) , 



(3.121) 



where n = ^. This leads one to 



TW^{su^'-Iju>') = -{V^u,)t>''' ~{TV^fi~E^)u>'-KfiE^Bi^^ j^^,) , (3.122) 



which is a typical starting point for the consideration of entropy production. 

We consider F^,^ as being first order in derivative, therefore the anomaly term in the 
basic hydrodynamic equation fl3.116p is of A^'th order in derivative. The left-hand side is 
V^j^, so it is expected that the first effects from anomaly should appear in at (A^ — l)'th 
order in derivative. This is a generalization of the case of A = 2 or in four dimensions 
where indeed the triangle anomaly affects j'^ at (A^ — 1) = 1st order. The corrections 
to the energy-momentum tensor are expected to be of higher order, presumably starting 
at A^'th order in derivative, which is beyond of our interest in this section. We stress 
that there should in general be many other non-anomalous viscous terms of (A — l)'th 
or less order in derivative, and we cannot possibly classify them all. However, the point 
is that the terms that are necessarily linear in the anomaly coefficients k do not mix with 
these other non-anomalous terms due to discrete symmetries (C, P) that we discuss in the 
previous section, so that it makes sense to discuss them separately. In 2A^ dimensions, 
the (C, P) charges are similar to those in four dimensions in the previous section, except 
55,^^,1) now has (C, P) = ((-l)^-\+l) and k has (C, P) = ((-1)^-^,-1), so that the 
transport coefficients are classified by there (C, P) charges as follows. 



To summarize, the transport coefficients of P = — 1 that are of interest for us are neces- 
sarily linear in k. 



(C,P) = (+1,+1) 
(C,P) = (-1,+1) 
(C,P) = ((-l)^-\-l) 
(CP) = ((-1)^,-1) 



f{T,p,,K) 



g{T, /i^/^^) 
fig{T, p^/^^) 

Kg{T, yU^K^) 

Kfig{T, p^/^^) 
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When one constructs the viscous corrections from various derivatives of thermody- 
namic quantities, the basic hydrodynamic equations fl3.116p can be used to remove some 
of the first order derivative terms up to higher order terms, so that one does not need 
them in constructing viscous terms at a given fixed order in derivative. As the total 
number of equations in f l3.116p is (A^ + 1), one expects to be able to remove {N + 1) first 
order derivative terms using the equations of motion. The general first order derivative 
terms are Vfj^u'^, F^^, and derivatives of any two independent thermodynamic scalars, say 
V^(T, or V^{p,n), because locally the plasma is completely specified by two indepen- 
dent thermodynamic scalars. It is a matter of choice which (A^ + 1) terms among the 
above first derivative terms are removed by using the equations of motion (13.116^ . For 
our convenience, we will remove 

V^p , u^V^^n , (3.123) 
so that the remaining available building blocks of constructing viscous terms are simply 

VX > F^^u > ^^"Vuu , (3.124) 

where A^'^ = u^u^ + g^^ is the projection operator to the space transverse to u^. 

To make things explicit, let us work out in detail how the above mentioned removal 
happens. Up to higher order derivative terms, the basic equations of motion f l3.116p are 
written as 

\/^f = M^V^n + nVX = ■ (3.125) 

From the second equation, one has 

uf'V = -n (y ^uf") . (3.126) 

On the other hand, contracting the first equation with Uy, one gets 

(e + p)(VX) = -«^V,.e , (3.127) 

and inserting this into the first equation of f l3.125p . one can obtain after an easy manip- 
ulation 

A'^^'Vf.p = nE"" - (e + p) u^'W ^.u'' . (3. 128) 
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Finally, from fl3.127p and writing e = e{p,n) and using fl3.126p . 

^'^'^ u^V^p-n(^\ (VX) , (3-129) 



dp I \dn 



so that one finally arrives at 



= L\ ■ (3-130) 



\dpJn 



The expressions ( I3.126p . f l3.128p . and fl3.130p indeed replace V^jo and u^V with V ^j,u" 
and F^j, up to higher order derivatives, so that we can remove them. 

We are interested in the anomaly- induced viscous corrections to the current v^. Since 
is a vector whose spatial component has P = —1 and we want a transport coefficient 
having P = — 1 to be linear in k, any anomalous viscous correction to should be a 
pseudo- vector whose spatial component has P = +1. It is easy to see that the first possible 
pseudo-vectors one can construct out of V ^u^ , F^,^, A^'^Vu'n and their derivatives indeed 
start to appear at (A^ — l)'th order in derivative, containing one e-tensor. It is also not 
difficult to see that the anomalous corrections to the energy-momentum tensor appear 
only at A^'th order and beyond. At (A^ — l)'th order of our interest, we find precisely A^ 
possible pseudo- vectors which can appear in and s'^j^_^y, 

%,t) = ^^'^'''^'---'''^'''''---^^^^^^ > s + t = N-l . 

(3.131) 

The B'^Qj^_^^ and q-, are 2 A^- dimensional generalizations of = ^e^'^°'^u^Fa/3 and 

(VaW/s) in four dimensions, while other -Bf^ ^ with < s < (A^ — 1) exist 



only in higher dimensions. One then introduces 2A^ transport coefficients ^(^s,t) and 



as 



^( A'' — 1 ), anomaly 



H ' ^(iV-1), anomaly = H ^(«,t)^&,t) ' (3.132) 

s+t=N-l s+t=N-l 



and inserts these into the entropy production formula ( ]3.122p to get some constraints on 
them. 



To proceed, we need the following formula derived from using the equations of motion 

{s + 1 
(e + p) 



= -77— Ti?&,) (V,P - nE,) - 2(Ar - 1 - .)i?,Pf,+i,_,) , (3.133) 
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where by definition B'^^ = if s > A^, and tlie above equality liolds true up to liiglier 
order corrections. 

Proof : Let us start from the definition 



V,i?f,,) = le^-^/^--^^^-- iy^u,) (V„,«^J ■ ■ ■ F,,,, ■ ■ ■ , (3.134) 



where other possible actions of V^j to (Va,W/3.) or F^.^. give simply zero using Bianchi 
identities of the Riemann tensor and the field strength tensor. It is most convenient to 
work in the local rest frame where = (z = 1,2,3) and V^ito = 0, so that the above 
becomes 



+ ^ '^ 6^-^^^-°^^-(V^^,)(V.,tx;,J---Fo,,--- . (3.135) 



Next, multiplying the equation of motion V^T^'' by 5;y(s,t) and using that u'^B^i^s t) = 0, 
one obtains up to higher derivitive corrections the following relation 

1 (e + p) ^M^.e'^A^,^,...,,^,... (V.,up,) ■ ■ ■ F,,s, ■■■ = Sf,,^) {-V,p + nE,) (3.136) 

which becomes in the local rest frame 



1 (e + p) e0-i/^--7i5i- (Vo^.,) (V„,i.^J ■ ■ ■ F,,s, ... = B^,^ {-V,p + nE,) . (3.137) 

The left-hand side is precisely proportional to the first term in (13.1351) . Finally, using 
-Po<5i = uqEs^ one manipulates the second term in (I3.135P as 

= ^o^,^e0^iM-i/3i-725.... (y^^^) (V„,M^J ■ ■ ■ F,,s, ■ ■ ■ (3.138) 
= -£;,^e^iOM-ift-7252-^^ (v^ii,) (V^.up,) ■ ■ ■ F^,s, ■■■ = -iVE^.Sf^Vi.^-i) • 

The fl37[35D . ^JSTh . and fl37[38|) prove our relation fl3:T33D . (QED) 

Using (13.1331) . the entropy production formula (I3.122p at (A^ — l)'th order upon in- 
serting (13.1321) becomes 

+ E fe(s,) + ,^i^'' -P(^.^) - 2T(^ - i^.i^r^,) ,(3.139) 
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where we formally define 

It is important to remember that the equations of motion relate E'^ to 

E^ = ^(A^'^V^p+(e + p)M'^VX) > (3-141) 

using (I3.128p . so one needs to be a bit careful when considering arbitrary possible inde- 
pendent configurations. The above entropy formula has a structure 

Km) + C(m)^m) ^(M) > (3-142) 

and the coefficients ^^(s,t) , involve only thermodynamic scalars and their derivatives. 
As -B(!jt) are clearly arbitrary and independent, the total coefficients in front of each of 
them should vanish in order to make sure positivity of entropy production, so the basic 
constraints one derives is in fact 

^M(M) + C(,,t)£^^ = . (3.143) 

An important observation is that E^ contains u'^ViyW^-piece while Af^i^s,t), C[s,t) are made of 
thermodynamic scalars and their derivatives only, without any u'^VyU^. Because u^VyU^ 
can vary independently of derivatives of thermodynamic scalars, imposing f l3.143p for 
arbitrary possible configurations gives us a simpler conclusion that 

^MM) = , q.,t) = , (3.144) 

separately. 

The reason why we expound on this subtlety is due to an interesting difference from 
the conformal case that we have studied in the previous section. Recall that in that case, 
E^ = ^"Dfj^p by equations of motion and the and C are also expressed in terms of 
of thermodynamic variables, so that one can no longer consider i?^ as being independent, 
and the constraint one has is only ( ]3.143p . not ( 13.144^ anymore. What happens in the 
conformal case is that u"V ^u^ and V^(T, Jj) are packaged to give only two independent 
combinations Vfj_{T,p,), not the general three, so that the number of constraints one gets 
is reduced. However, what saves us is the conformal symmetry constraint on the transport 
coefficients / = f{p,) which is given from the start, so that one is still able to solve the 
system of reduced constraints, as we do in the previous section. 
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Therefore, one has the system of 2N algebro-differential equations 

( V^/i) ecM) - + ( V^p) A.,) = , (3.145) 

( s -L i)Tn 

^is,t) + -j^^^D^s,t)-2T{N-s)D^,.,,t+,^ = , (3.146) 

with the boundary condition D(_i^7v) = i^'^A*- We can solve them analytically. Guided 
by Ref . [33] , one chooses (p, p) as the two independent thermodynamic scalars describing 
the plasma, so that any transport coefficient is considered as a function of {p,p.)- Note 
that this is only a convenient choice of parameters. From fl3.145p . we then have 

and using the fact 

/<9T\ T 



p 



dp J. (e + p) 



(3.148) 



which can be derived by considering 



_ (e + p) 

dp = sdT + ndfj, = sdT + nd (Tfi) = {s + n^) dT + nTdfi = — — — dT + nTdfi , (3.149) 
the first equation gives us 

which is solved by 

D^s,t)=T'^'D^s,t){f^) . (3.151) 
Then, j) is given by the second equation of (13.1471) and 

fdT\ nT^ 



d^il„ (e + p) 



(3.152) 



which can be obtained in a similar way. The remaining step is to find -D(s,t). Inserting 
(13.1511) into (13.1461) and using some of the above equations, one arrives at 



-is+1) (^^^ -2{N-s) T^+iD(,_i,t+i) 

-2{N-s)T'+'D^s-i,t+i) = , (3.153) 



-1 dD^ 
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which simply results in an iteration equation 

with an initial condition i)(_i^Ar) = j^nft. This is easily solved as§ 



A-.(f) = (7TliiF:7b^''^'« . (3.155) 



From (13.1511) and the second equation of (13.1471) . one finally has 

D = \Y\ , (3.156) 



3.2 AdS/CFT correspondence 

We will close this section by confirming our results in the previous subsection in the 
AdS/CFT correspondence via fiuid/gravity computation |59]. We start from the holo- 
graphic bulk action in [d + 1) = {2N + 1) dimensions, 

167rG(,+i)) C = R + did-l)- ^Fm^F*^^ - -^eA-^^iQi-^-Q-^^^^^^^^ . . . p^^^^ ^ 

(3.158) 

whose equations of motion are 



1 " FPQ] 1 



Rmn + yd+ ^ ^^ FpqF j qmn — -FMpFjy — , 
9 



The static charged black hole solution is 

d-l 2 

ds"^ = -r^V{r)df + 2dtdr + ^ (rfx*) , A = — j-^dt , (3.160) 



1=1 

with 



nr) = l-^+ . (3.161) 



^We again neglect possible integration constants which are not fixed by the method. 
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For fluid/gravity computation, the general boosted solution is 

(iS(o) = —r'^V{r)u^Uudx^dx''—2Ufj,dx^dr+r'^ [g^i, + u^Uu) dx^dx" , ^4(0) = -^^u^dx^ . 

^ (3.162) 

and one lets the parameters of the solution (m(x), q{x), u^{x)) to vary in space-time, and 
systematically adds corrections to the above zero'th order metric in derivarive expansions 
of {m{x) , q{x) , u^^ {x)) in order to satisfy the equations of motion f l3.159p . For our case, 
we also need to introduce the external gauge potential in the derivative expansion, so we 
in fact should extend the zero'th order bulk gauge field as 

Ao) = ^M^)dx^ + M^)dx^ ' (3-163) 

where yl^(x), by an abuse of notation, is understood as an external gauge potential. The 
Ansatz is then 

ds^ = dsfo) + g^M'^dx^^dx^ , A = A^o) + A^^dx^^ , (3.164) 

k 

where the corrections of fc'th order ( (7^^^^ , A^'' ) contain total k number of derivatives of 
{'m{x),q{x),u^{x), A^{x)). It is convenient to use general coordinate/[/(l) gauge trans- 
formations to work in the gauge such that 

9ir^ = , gi'J-^u, , E^?J? = , 4'^ = . (3.165) 

i 

At order k, the equations for {g^^\, A^j^) are simple second order linear ODE's along r 
since x'^ derivatives of them are of higher order, and the sources for these linear ODE are 
given in terms of solutions up to {k — l)'th order and their derivatives, and of order k in 
total number of derivatives. It is an important fact that the linear ODE operators acting 
on {gM^N, A^j^) are universal, independent of k and can thus be determined at the first 
order. See Ref.[59] for details. 

To find equations for {g^j^^j^, A^^j) at the position, say = 0, it is most convenient to 
work in the local rest frame where u^^ = (—1,0) at x'^ = 0, and one can locally expand 
{'m{x),q{x),u^{x), A^{x)) in powers of x'^'s up to the order of interest. For example, at 
first order in derivatives one can use 

u^{x) = {-l^x'^df.Ui) , {m{x),q{x)) = {m,q)+x'' (d^m,df,q) , A^{x) = -^F^^x" 

(3.166) 
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One can classify {q'mni ) according to their representations under the spatial symmetry 
group SO{d — 1) of our local rest frame, and since we are interested in the dissipative 
transverse U{1) current, only the vector modes {gii\Af'^) are relevant. 

It is easy to check that the Chern-Simons term in the equations of motion fl3.159p 
indeed starts appearing in the procedure at order k = (N—l) in the vector mode equations 
due to the total anti-symmetrization of e-tensor and the fact that the only non-zero zero'th 
order -F^/jv is F^^^ . The linear differential operators along r acting on {gl'^\A\^^) can be 
found easily from Ref . [36] . Since we are aiming to effects from anomaly or Chern-Simons 
term, we keep only source terms that come from the Chern-Simons term. The relevant 
equations to solve are 

a. f/+^aJ^))+(d-2)g(a.,Ar"'0 = ' (3.167) 



(7V-l)^ 



a,(/-V(r)(a,Af^'))) + (d-2)ga, l^^j = S^ir) , (3.168 



where the source Si{r) coming from the Chern-Simons term is found to be 

S^{r) = -^d {d' - 4) Cq^^^'^^'--^--^'--^F^Xi^) ■ ■ ■ F^'l^.^Jr) , (3.169) 
where 

F^l\r) = ^ id,u, - d,uj) + F,j . (3.170) 

Recalling = — 1 in our rest frame and expanding products of Fjl'^ in the source Si{r), 
one easily recognizes ^^-structure appearing. Explicitly, one has 

^«(^) = -^)CgT. iH' J),.-, ^iM) > (3-171) 

It is not difficult to solve fl3.167p and fl3.168p [36] ; first integrate (13.1680 to get 



(3.172) 

where considering the boundary condition at the horizon r = rn, the integration constant 
Ci is equal to gti^~^\rH) which should be determined later. Inserting this into (I3.167P 
removing A^^ ^\ one obtains a second order differential equation for g[f which turns 
out to be after some manipulations 

dr\r'^'{V{r)fd,{^,\ \ =I{r) , (3.173) 
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which is readily integrated as 



ft'f-'V) = rV(r) y^d/^-^^j-^-^ [l/r"'(r") - r%-V(r„)Q j . (3.174) 



where integration constants are fixed by considering regularity boundary conditions at 
the horizon. Then, A-^"^^ is given by integrating fl3.172p . 

The above is the complete solution for {g['^'~^\ Af^~^^) except Q still needs to be fixed. 
The solution is regular for any Ci and what fixes it is the Landau frame condition that 
Tti^ = in our local rest frame, or equivalently ^:a^-piece in the near boundary asymp- 
totics of g[f should vanish. This brings us the condition for Q, 

/ drI{r) = r''j^^V\rH)C, , (3.176) 



which finally determines Ci as 



{d-2)qr% foo 1 



Ci = ---r^— dr-- dr'S,{r') . (3.177) 



This completes the solution for [g^^ ^^)- 

Near the boundary r — > oo, a[^~^^ has the asymptotics. 



and from this the U{1) current from anomaly at (A^ — l)'th order is obtained as 

_ {d-2) /-2^(^-i). ^ _ 1 ( rdrS (r) + ^±L^C- 



'(d+l) '^^ iU7rLT(d+i) \Jr„ rjj- 



(3.179) 

There is one point we need to make clear; in fact, the full boundary current gets an 
additional contribution from the variation of the bulk Chern-Simons term as was first 
stressed by Ref . [12] . 



(Ar_l) _ (Ar_i) Cd rf,ui5i-,y,^_-LSr,.i A p ^ ...p ^ ("^ T SOl 



,(N-l) _ ,,{Af-l) ^d 

r{d+l) 

However, rigorous holographic renormalization of the theory in Ref. [00] shows that what 
enters in the energy-momentum Ward identity 

Q^T^u^p^u^^ , (3.181) 
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is the current obtained only by the near boundary asymptotics hke z^l^^nomaiy ^^e 
above, not the full current ju^mi including additional contribution from the Chern-Simons 
term The basic reason for this is that the Chern-Simons term is topological and does 
not couple to the metric at all, so that it gives no contribution to the energy-momentum 
tensor. The choice of a current one is working with does not matter as long as one is clear 
about its definition and the correct Ward identities. Since our discussions are based on 
the energy-momentum Ward identity of the form f l3.18ip . the correct current we have to 

use here is Z/^,anomaly 

Finally, we need to fix the coefficient C in front of the Chern-Simons term, and the 
point in the above paragraph is also relevant here. Recall that our basic hydrodynamic 
equations (13.1161) assume the energy- momentum Ward identity of the form f l3.18ip . so 
that we need to match C such that the divergence of the current obtained only from the 
near boundary asymptotics, not of the full current, agrees with the second equation of 
(13.1161) . The bulk equation of motion of the ^7(1) gauge field (I3.159P near the boundary 
r — 7- oo dictates that 

{d + 2)C n^ui--fiNUN p . . .p fQ 1 fiO^ 

327rG(d+i) 

where j'^ in the above is the current obtained only from the near boundary asymptotics 
as in (I3.179p . Comparing with (I3.116p . we can fix C as 

d\d + 2)^ ■ ^^-^^^^ 
With f lXTTT]) . f lXTTT]) . flXTT^ . and flCT^ . one finally obtains after some algebra 



(iV-l) 

i, anomaly / ^ ^(s^t)-'-' {^s,t) ' 



Ee(.,*)4,) , (3.184) 



where the expressions for (^(^^j) is 

F 2-(iV-l)! ^g-+i / (s + l\(d-2\ 

(s + l)!(iV-s-l)!^g+i)('^-2) 1^^ \s^2)\ d '^^'^-2 1 • ^^-^^^^ 



mr 



H 



Using the relations 



, dm {d-2)q q (n^Qa\ 

e + p = dp=—— , n = — — , ^=^Z2 , (3.186) 

167rG(rf+i) 167rG(d+i) 



H 



^We thank A.Yarom for discussions on this. 
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it is easily checked that this agrees precisely with our result fl3.157p in the previous sub- 
section. 
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